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ABSTRACT 

Ultra- high energy cosmic rays (UHECRs) accelerated in the jets of active galactic 
nuclei can accumulate in high magnetic field, ~ 100 kpc-scale regions surrounding pow- 
erful radio galaxies. Photohadronic processes involving UHECRs and photons of the 
extragalactic background light make ultra-relativistic electrons and positrons that ini- 
tiate electromagnetic cascades, leading to the production of a 7-ray synchrotron halo. 
We calculate the halo emission in the case of Cygnus A and show that it should be 
detectable with the Fermi Gamma ray Space Telescope and possibly detectable with 
ground-based 7-ray telescopes if radio galaxies are the sources of UHECRs. 

Subject headings: Galaxies: individual (Cygnus A) — cosmic rays: theory — galaxies: 
active — galaxies: jets — radiation processes: nonthermal 



Introduction 



Active galactic nuclei (AGN) and gamma-ray bursts are considered as two of the most plausible 
class es of astrophysical acce lerators of extragalactic ultr a- high energy cosmic rays (UHEC Rs; see, 



e.g 



Halzen Sz Hooperl 12002 ) . The recent report of the iPierre Auger Collaboration! (|2007l ) about 
clustering of the arrival directions of UHECRs with energies £>6x 10 19 eV within ks 3° of the 
directions to AGN at distances d < 75 Mpc strongly suggests that effective production of cosmic 
rays with energies E ~ 10 20 eV takes place in at least one of these source classes. Because of 
photohadronic GZK interactions of protons or ions with the CMB radiation, the study of super- 
GZK UHECRs from sourc es at d > 100 Mpc becomes impossible with cosmic-ray detectors like the 
Pierre Auger Observatory (jHarari et al.l l2006). Powerful AGN, as well as GRBs, are mostly located 
at larger distances. 

Relativistic beams of energy from the central nuclei of AGN are t hought to powe r the multi-kpc 
scale radio lobes of powerful galaxies and form an extended cavity (jScheuerl 11974 ) . Acceleration 
of UHECRs in the compact inner jets of the radio galaxy on the sub-parsec scale, followed by 
production of collimated beams of ultra-high energy neutrons and gamma-rays, provides a specific 
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mechanism to transport energy to the radio lobes and cavity (jAtovan &: Dermerll2003l ). When the 
neutron-decay UHECR protons interact with the extragalactic background light (EBL), which is 
dominated by the CMB radiation, ultra-relativistic electrons (including positrons) and 7 rays with 
E > 10 18 eV are produced. Such secondaries can initiate pair-photon cascades to form large multi- 



Mpc scale halos of GeV /TeV radiation du e to Compton and synchrotron processes (jAharonian et al 



1994; 


Aharonian 


2002; 


Inoue et al. 


2005) 



In this Letter, we predict that synchrotron GeV fluxes from UHECR AGN sources are de- 
tectable with the Fermi Gamma ray Space Telescope (FGST; formerly the Gamma ray Large Area 
Space Telescope, GLAST) if UHECRs are captured in the vicinity of radio galaxies for sufficiently 
long times. Magnetic fields at the > fiG level in the kpc - Mpc vicinity from the AGN core are 
required to isotropize UHECRs accelerated by jets of radio galaxies. Indeed, for protons with en- 
ergy E = W 20 E 2 o eV, the mean magnetic field B required to provide gyroradii smaller than size r 
is B > 10~ 4 £ , 20?"kpc G> where r^ pc = r/lkpc is the spatial scale where the isotropization occurs. 

Here we consider the specific case of the powerful radio galaxy Cygnus A, where the mean 
magnetic field in the surrounding cavity could reach 10 - 100 at 100 kpc scales. Its properties 
are considered in Section 2, and calculations are presented in Section 3. We summarize in Section 
4. 



Model Parameters 



Magnetic fields B > 1 uG could be prese nt at < 1 Mpc scales in clusters of galaxies (IKim et al 



1990 



Feretti et al 



1995 



Ferrari et al.ll2008l ). Even higher magnetic fields, B ~ 10 4 G, could be 



present in the r < 1 00 kpc vicinity of c D galaxies near the center of galaxy clusters like the powerful 
radio galaxy Cyg A (jWilson et al.ll2006l ). In such magnetic fields, synchrotron radiation of electrons 
with E ~ 10 18 eV produced in photomeson interactions is in the TeV domain. Furthermore, 
synchrotron radiation of lower energy electrons produced in p + 7 — > p + e + + e~ interactions by 
UHECRs extends to the GeV domain and could become detectable with the Fermi Telescope. Here 
we consider the detectability of these fluxes from the powerful and well-studied radio galaxy Cyg 
A. 



2.1. Power of the Cosmic Ray Accelerator 



Observations of Cygnus A (z = 0.056, luminosity distance ~ 240 Mpc) by the Chandra X-ray 
Observatory show that the central ~ 60 kpc x 120 kpc size luminous prolate spheroidal region of 
Cygnus A, called a "cavity," can be understood as a shock expanding into the accretion cooling- 
flow gas (Wilso net al.ll2006l ). The kinetic power of expansion of this X-ray cavity deduced by 
Wilson et al.l ( 20061 ) is L exp ~ 1.2 x 10 erg s . This is much lar ger than the total r adio luminosity 



L 



lobes 



10 45 erg s 1 of the two bright radio lobes of Cygnus A ( Perlev et al.lll984l ). The value of 
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L exp should be considered as a lower limit to the overall power Lqr of cosmic rays injected into 
the cavity, if the cosmic-ray power is assumed to drive the expansion of the cavity. A few times 
larger power, Lqr ~ 4 x 10 46 ergs s _1 therefore seems a reasonable assumption. 



This value is in a good agreement with the "neutral beam" model (jAtoyan &: Dermerll200ll . 



2003), which explains the collimated relativistic X-ray jets that remain s traight up to ~ 1 Mpc 
distances in sources such as Pictor A (jWilson. Young, and Shopbeilll200ll ) as the result of energy 
transport by beams of UHE neutrons and gamma-rays. These linear X-ray features, surrounded 
by a broader and less collimated radio structure in P ictor A, and coincident with narrow radio 
struc tures exhibiting bends and deflections in Cyg A (jCarilli et al.lll996l ; ISteenbrugge fc Blundell 
20081 ). terminate in X-ray hot spots. Because of the large inclination angle, the X-ray jets in Cyg A 
cannot be detected. However detection of bright X-ray hot spots located at ~ 60kpc distances on 
opposite sides of the nucleus (jWilson et al.ll2006l ) strongly implies production of collimated X-ray 
jets also in Cyg A. 

Detailed calculations in the framework of this model show that neutral beams of ultrarelativistic 
neutrons and 7 rays, produced by the compact relativistic jets in the central sub-parsec scale 
environment of FRII galaxies, can take a few percent of the total power of this inner jet. That 
energy is then deposited, after /3-decay of neutrons n — > p + e + v e , on distance scales Id ~ E20 Mpc, 
and also via photopair production of UHE gamma-rays through the process 7 + 7 — > e~ + e + . 
These secondary charged relativistic particles initially form a beam in the same direction as the jet, 
and can effectively interact with and transfer momentum and energy to the ambient magnetized 
medium, and thus can be the basic energy reservoir for the jet at large dist ances from Cyg A. The 



total radio lumino sity ~ 9 x 10 44 ergs 1 detected from the lobes of Cyg A ( Carilli &: Barthel 



1996 



Wilson et al.ll2006l ) imposes the minimum power requirement for the beam. Given the ~ (2 — 3)% 
efficiency of neutral beam production, the acceleration power of the CRs in the relativistic inner 
jet must be L C r > (3 - 5) X 10 46 erg/s. This power is then released in CRs when the inner 
jet decelerates to subrelativistic speeds in the dense medium at kpc scales. This scenario is in 
agreement with th e assumption that the expansion of the cavity against the cooling flow observed 
(jSmith et al.ll2002l ) at distances < 70 - 100 kpc is powered by cosmic-ray pressure. 

Protons are accelerated in the inner jet of Cyg A to a maximum energy E max ~ 10 20 eV, 
consistent with size scales and magnetic fields inferred from a synchrotron model of blazars corre- 
sponding to Cyg A if observed along its jet. This value of E max is in accord with the spatial extent 
of the X-ray hot spot where the large-scale jet terminates due to the /3-decay of the remaining high- 
energy neutrons in the beam. At this distance, the injection power from the decaying neutrons is 
balanced by the ram pressure of the external medium. Thus, for Id ~ 100 kpc-long jets originating 
from neutral beams, the inner jets must accelerate protons to E > 10 19 eV. In our model, we assume 
that acceleration of UHECRs to E ~ 10 20 eV occurs. 
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2.2. Magnetic Field Strength 



The lower limit to the equ ipartition magneti c field in the radio lobes of Cyg A inferred from 



radio observations is « 60 /iG (jCarilli et al.l 119911). The equipartitio n magnetic field in the cavity 



derived from X-ray observations is B ~ 120 /iG (jWilson et al.l 120061 ) . Moreover, the field outside 
the cavity could also be very high given the large thermal electron energy density inferred from 
the X-ray emission at < 8' or in the < 500 kpc environment around Cyg A. The pressure of hot 
thermal gas decreases from 2.3 x 10 



io er g cm 3 j n regions near the c avity to < 10 11 erg 



cm -3 at distances > 500 kpc from the center of Cyg A ( Smith et al.ll2002l ). The corresponding 
equipartition magnetic fields would then vary from 
region. 



80 fiG to ~ 15 /xG at the periphery of this 



2.3. Injection Age 



The duration of injec tion of UHECRs, i.e ., the jet injection age, represents one of the important 
parameters of the model. I Wilson et al.l (|2006l ) find that the expansion age of the cavity, determined 
from the speed of the shock deduced from the analysis of X-ray data, is t ~ 3 x 10 7 yr. The injection 
age can be larger than this dynamical age because the derived value neglects the magnetic-field 



pressu re of the intracluster medium (ICM) upstream of the shock, as recognized by iWilson et al 
(|2006l ). Also, this age estimate neglects the infall velocity of the cooling flo w. For the total 
mass M ~ 2 x 10 13 M Q enclosed at r < 50 kpc distances ( Smith et al.1 120021 ) . the virial speed 
Vvh ~ c/^/r/r s 9* 2000 km/s, where r s = 2GM/c 2 
radius for the total mass M. The value of f v ir 



6 x 10 18 cm is, formally, the Schwarzschild 
gives a measure of the accretion/cooling flow 



vel ocities at rad i us r, and is comparable to the average ~ 1500 km/s speed of the shock derived 



by IWilson et al.l (|2006l ) in the rest frame of the fluid upstream of the shock. These factors can 
significantly decrease the speed of expansion of the cavity in the stationary fr ame, so that t he rea l 
injection age of the cavity can be significantly larger than the age inferred by I Wilson et aJj ([2006). 



We now estimate the age of activity of the black-hole jet from energetics arguments. The 
inferred jet power, Lj e t — 4 x 10 46 erg/s, correspo nds to 12% of the E ddington luminosity for a 
black-hole mass Mbh — 2.5 x 1O 9 M0 in Cygnus A ( Tadhunter et al.ll2003l ). To produce such power, 
the black hole should accrete mass at the rate M = Lj^/rjc 2 = (7 /t7_i)Mq yr -1 with an efficiency 
i] = 10 x ?7 i , with ?7„ i ~ 1. The age of the central black hole is estimated by its growth time 
t ~ Mbh/M = 3.6 x 10 8 ?7_i yr, giving an upper estimate for the jet's age as the jet might be active 
for only a fraction of the BH growth phase. These estimates are in accord wit h the characteristic 
j et ag e t ~ 10 8 yr inferred from a model for "cocoon" (or cavity) dynamics by lBegelman &; Cioffi 
dl989l). 
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Fig. 1. — Comparison of timescales for energy losses through photomeson (dashed curve) and 
photopair (dot-dashed curve) production, confinement (solid line) with B = 20 fiG and r = 500 
kpc, and injection with t age = 10 8 yrs (triple-dot-dashed line). The inset shows the fraction A(E) 
of energy extracted by photohadronic processes from UHECRs with energy E. 



2.4. Cosmic Ray Diffusive Confinement Time 

The maximum confinement timescale t con f of UHECRs in a source of size r is given in the Bohm 
diffusion approximation by t con f = r 2 /2L>B> where the Bohm diffusion coefficient Db = and 
ri = E/QB is the Larmor radius of a particle with charge Q. From this expression we obtain the 
UHECR proton confinement time 

t conf ^ 0.95 x 10 7 (r/100 kpc) 2 ^ 1 (5/20 fiG) yr . (1) 

For r ~ 50 kpc and B ^ 120 fiG, t conf 15 Myr for « 10 20 eV protons. 

Fig. [H compares timescales ti oss for energy losses due to photomeson and photopair production 
with values of t con f for characteristic parameters r = 500 kpc and B = 20 fj,G in the region surround- 
ing Cyg A. Also shown is the value of t age = 10 8 yrs for the assumed injection age. The fraction of 
energy of accelerated protons that could be extracted is given by A(E) = mm(t con f,ti oss , t age )/ti oss . 
This figure shows that > 3% of the total energy of UHECRs with i?>3x 10 18 eV can be extracted 
through photohadronic processes. Because of the increased confinement time and the much lower 
production threshold for photopair than photomeson processes, the photopair process can make 
a comparable or dominant contribution to the electromagnetic channel compared to photomeson 
processes. 
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Fig. 2. — Fluxes from the electromagnetic cascade initiated within 0.5 Mpc distances from the core 
of Cyg A, assuming an average magnetic field B = 20 /uG in the cavity. The total injection power 
of UHECR protons is 4 x 10 46 erg/s. The photon fluxes expected for t age = 30Myr, 100 Myr and 
300 Myr are shown by the open dots, full dots and stars, respectively. The solid, dashed, and dot- 
dashed curves show three generations of synchrotron (heavy curves) and Compton (light curves) 
cascade radiations in the case of t = 100 Myr. Note that the third generation of Compton cascade 
radiation is too weak to appear on the plot. Sensitivities for the FGST and VERITAS are also 
shown. 



3. Results 

Fig. [2] shows 7-ray fluxes calculated for the electromagnetic cascade initiated by the injection 
of 4 x 10 46 erg s _1 of UHECR protons into the cavity, and further cosmic ray interactions with the 
EBL in the r = 0.5 Mpc region surrounding Cyg A. Even though the EBL is dominated by the 
CMBR, interacti ons with the dif f use in frared /optical radiation are also included in the calculations, 
using the EBL of iPrimack et al.l (|2005l ). The injection spectrum of UHECR protons is a power-law 
with an a = —2.1 number spectral index with a low-energy cutoff at E = 1 TeV and a high-energy 
exponential cutoff at E = 3 x 10 20 eV. We assume a mean magnetic field B = 20 fiG. Escape 
of particles is given by the Bohm diffusion a pproximation in a s i ngle z one model. The method of 
calculation follows the approach described by lAtoyan &: Dermerl (120031 ). 



The solid, dashed, and dot-dashed curves in Fig. [2] show the contributions to the fluxes from 
the photohadronic secondary electrons and the first two generations of cascade electrons for in- 
jection ages t age = 10 8 yr. The secondary electrons also include the electrons from 7r°-decay 7 
rays that convert promptly into electron-positron pairs inside the source. In the strong magnetic 
field of the Cyg A environment, the Compton flux is dominated by synchrotron radiation from the 
photohadronic secondaries and the first generation of cascade electrons. 



The open dots, full dots and stars in Fig. [2]show the received E-F(E) spectral energy fluxes for 
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injection ages t age = 30 Myr, 100 Myr and 300 Myr, respectively. The lower and higher energy peaks 
in the spectral energy distributions primarily result from photopair and photomeson processes, 
respectively. The sensitivities for a one-year, 5a detection of a point source with the FGST in the 
scanning mode0 and for a 50 hour, 5a detection with VERITAS are shown. The lower bound of 
the VERITAS sensitivity applies to a point source, and the upper bound to a source of 25' extent. 



4. Discussion and Conclusions 

If the radio lobes of Cyg A are powered by UHECR production from the inner pc-scale jets, then 
trapping of these particles in the surrounding strong magnetic-field region leads to the production 
of secondary 7 rays that should be significantly detected with the FGST in one year of observation 
if the jet injection age is > 100 Myr. If radio galaxies are not the sources of UHECRs, then Cygnus 
A will not be detected by the FGST. Cyg A could also be detected with VERITAS in a 50 hour 
pointing, depending on the duration of activity of the central engine and the level of the EBL. 

Detection of GeV 7 rays from Cyg A with the FGST might also be expected to arise from 
other processes. The ~ 10 - 100 GeV radio-emitting electrons fr om the lobes of radio galaxies 



will Compton-scatter CMB photons to MeV - GeV energies (e.g., ICheungj 120071 ). For the strong 
magnetic field, ~ 60 fiG, in the lobes of Cyg A, however, the ratio of the magnetic-field to CMBR 
energy densities is ~ 400. Thus the total energy flux of Compton-scattered CMBR from Cyg A 
is 10 45 erg s" 1 /[400(4vr(i 2 )] ^ 4 x 10~ 13 ergs cnr 2 s -1 , with the E ■ F(E) flux a factor of ~ 5 
- 10 lower. As can be seen from Fig. (2j this process is almost two orders of magnitude below the 
UHECR-induced synchrotron flux, and well below the FGST sensitivity. 



Inoue et al.l (|2005l ) considered fluxes expected from the < 1 Mpc halos of clusters of galaxies 
with weaker magnetic fields, B ~ 0.1-1 [i G in a model where acceleration of UHECRs occurs 
in accretion shocks in the cluster. Because of lower maximum energies of accelerated protons, 
E < 10 19 eV, and lower m agnetic fields, this model pre dicts hard spectral fluxes of Compton origin 



peaking at TeV energies. iGabici &; Aharonianl (|2005l ) predicted that synchrotron radiation from 
> 10 18 eV electrons is produced by secondaries of UHECRs that leave the acceleration region and 
travel nearly rectilinearly through weak intergalactic mag netic fields at the level B ~ 10~ 7 - 10" 9 G. 
These sources would appear as point-like quiescent GeV - TeV sources with spectra in the GeV 
domain as hard as a = —1.5, and very soft, a < —3 spectra in the 100 GeV - TeV domain. 

In contrast to both these models, we predict soft 0.1 - 1 GeV spectra with a = —2.5 and hard, 
a = —2 spectra at TeV energies due to the much higher magnetic field in the confinement region. 
These models can be clearly distinguished if Cygnus A is resolved by the Fermi Gamma ray Space 
Telescope or the ground-based 7-ray telescopes, as the emission region in our model subtends an 



www-glast.slac.stanford.edu/software/IS/glast_lat_performance.htm; note that the FGST sensitivity will be some- 
what poorer than shown at the 5.76° galactic latitude of Cyg A due to Galactic background. 
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angle « 10'. 



Because Cygnus A lies outside the GZK horizon, only UHECRs with energy below the GZK 
energy could be correlated with this source. Other closer FRII radio galaxies that are correlated 
with the arrival directions of UHECRs are, however, candidate sources of 7 rays made through the 
mechanism p roposed here. IGR J21247+5058 at z = 0.02 or d ~ 80 Mpc, recently discovered with 
INTEGRAL (jMolina et alj|20071 ). is 2.1 degrees away from a HiRes Stereo event with E > 56 EeV 
(C. C. Cheung, private communication, 2008)1 PKS 2158-380 at k, 1 40 Mpc is also within 3.2° 
degrees of an Auger UHECR with E > 57 EeV (jMoskalenko et al.ll2008l ). By comparison with Cyg 
A, these are low luminosity FRIIs, and their predicted flux level will require detailed modeling for 
each source, as done here for Cyg A. Variability of 7-ray emission would rule out our model. 
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